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Abstract 

Fermi Ball is a kind of nontopological soliton with fermions trapped in its domain wall, and is 
suggested to arises from the spontaneous symmetry breaking of the approximate Z2 symmetry in 
the early universe. We find that the neutral thin-wall Fermi Ball is stable in the limited region of 
the scalar self-coupling constant A and the Yukawa coupling constant G. We find that the Fermi 
Ball is stabilized due to the curvature effect of the domain wall caused by the fermion sector. We 
also discuss whether such stable Fermi Balls may contribute to the cold dark matter. 

PACS numbers: 05.45.Yv, 95.35. +d 
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I. INTRODUCTION 



In quantum field theory and especially in the cosmological context, various models have 
been discussed where the spontaneous breaking of a discrete symmetry produces domain 
walls fl], 3 . If the symmetry is approximate and broken spontaneously in the history of the 
expanding universe, the false vacuum regions shrink due to the energy density difference 
0] and this process is accelerated by the surface tension. The shrinking will stop before 
collapsing if there exist zero mode solutions for the fermions bound in such domain walls 
0, Q and the Fermi pressure of such fermions becomes comparable to the shrinking force 
due to the surface tension and the volume energy. These objects, a kind of non-topological 
solitons, are called Fermi Balls and introduced as a candidate for cold dark matter |P, ^ . 
They were also suggested in the baryon-separation scenario with the QCD energy scale [§. 
However, the stability of the Fermi Balls which should be essential to be a candidate for 
cold dark matter has not been fully examined. 

We consider in this paper the stability of the thin-wall Fermi Balls against the fragmen- 
tation. The Fermi Ball, if it is electrically neutral as was firstly proposed |p, is unstable 
against the fragmentation, when the domain wall thickness 6b can be neglected in compar- 
ison with the Fermi Ball radius R. In the previous paper [Q, we estimated the corrections 
caused by the finite 6fy effect (corresponding to the curvature effect of the domain wall), and 
found that the neutral Fermi Ball is unstable even if the curvature effect is included. Though 



the electrically charged Fermi Ball introduced by Morris |]T0| has the stability against the 
deformation from the spherical shape |rU|, they are metastable ^ and not absolutely 
stable against the fragmentation. We do not discuss the charged Fermi Ball in the present 
paper. 

In the above discussions, though, the thickness 6f of the distribution of the fermions 
confined in the wall is much less than the wall thickness 6^ and can be neglected. This 
approximation is valid only for the case where the fermions are tightly bound to the domain 
wall. We in this paper consider the case where 6f cannot be neglected. In order to esti- 
mate the corrections caused by the finite 6f effect, we use the perturbation approximation 
expanding the scalar field 0, the fermion field ip and the static energy of the Fermi Ball E 
in the power of 6b/ R. We examine the stability of the Fermi Ball against the fragmentation 
at each level of the perturbation. 
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The contents of the present paper are organized as follows: We first explain the neutral 
Fermi Ball model to clarify terminologies and method of energy estimation in Sec.|I|. We 
next examine the stability of the Fermi Ball within the leading (zeroth) order of Sb/R- 
perturbation but with finite 5/ in Sec.pTT|, and regain the result obtained in the previous 
works in the limit 5j ^ ||^. The effect of the finite Sj in the higher order corrections 
of 5;,/-R-perturbation is investigated in Sec.|V[ We discuss constraints on the Fermi Ball 
parameters from the experimental viewpoint in Sec.0. We summarize the obtained results 



in Sec. VI 



II. METHOD FOR ESTIMATING FERMI BALL ENERGY 

We consider the simple model with the Lagrangian density, 

jC = ^id,^f + ^itrd,-G<P)^-Ui<P) , (1) 

where the scalar potential U{(f)) is approximately double-well shaped, 

W) = ^(0'-^')' + A(0) . (2) 

Here, the first term has the Z2 symmetry under (f) — </>, and the second term violates the 
symmetry though it is assumed to be much smaller than the first one, A ~ \A{v) — A{—v) \ ^ 
Af^. Supposing A(— f) > A(f), we call the region with cj) = v and that with cf) = —v as the 
true vacuum and the false vacuum, respectively. In the following, we neglect the second term 
except for the explicit discussion in Sec.^. Fermi Ball is the ground state of the system 
with the total number of the fermion being fixed: 

Nf = J d^x . (3) 

The classical fields (l){x,t) and \l/(a;, t) for Fermi Ball extremize 

L[(j),-^;ef] = Jd^x C + ef(^Jd^x^^-^-Nf^ , (4) 

with ef the Lagrange multiplier. The static fields thus satisfy 

{ap + G4>p)^ = ef^ , (5) 
-V20 + + ^0(02 - = , (6) 



where a = 7^7, /3 = 7° and p = —iV. 

Assuming the spherical symmetry of we take it as a function of the radial coordinate 
r. We first consider Eq.@. Let "^{x) be the eigenfunction of J , Jz and P: 

Then, is written as 



(8) 



where yij{0, if) and yj/jid, are the spherical spinors having the eigenvalues J = l — uj/2 = 
r + uj/2 and M, with u = ±1. We take 3^;^5 = i^^/r)yu and note P = (-1)' = {-iy+^^^. 
Substituting Eq.(|D into Eq.(|D, we get 



(9) 



{pr + tf)f-G<Pg = ejg 

where Pr = —i{d/dr) and K = uj{J + 1/2). The equation is compactly written in terms 
ofV^(r)^g): 

Hf^ = ef^ , (10) 

with 

K 

Hf = (JiPr + (T2— + (TsGcj) . (11) 

The radial coordinate r is hereafter replaced hj w = r — R. 

We next consider Eq.(D. Noting = Ekm^m'^m = i^/4^T^r'^)Y.Ki'^\K\)^^a3ip, we 
have 

^ + = ^0(0^ - v^) + T71^2 E(2|^I)^V3^ . (12) 

dw^ R + w dw 2 4:n{R + ty)^ 



The energy of the Fermi Ball is expressed in terms of and ip as follows: 

E = Ef + Eb , (13) 

where Ef is the fermi energy, 

Ef = J d^x {ap + G(j)(3) ^ 



/+00 
dw^^Hftfj , (14) 



and Eh is the surface energy |23], 



hoo 



oo 



Note that we estimate the energy by integrating the above integrands not from —R but from 
— oo, since most of the contribution comes from the region near the surface and the error 
due to this approximation is exponentially small (oc ^-Const.vRy 

III. LEADING ORDER OF t^t/i^-PERTURBATION 

Let us expand the fields, 



= 00 + 01 + • 



(16) 



and the Hamiltonian, 



Hf = Ho + H, + H2 + --- , (17) 

K 

'r 

■ „ '" 

in the power of 5h/R (we need the expressions (|13|) to ([T5|) in the power of 5h/R up to 



K 

Hq = aiPr + 0-2— + (T3G0O , (18) 

K 

Hi = -(T2^w + asG^i , (19) 
H2 = a2^w^ , (20) 



the leading and the next-to- leading order). We also expand E = J2iEi, Ef = J2iE^f\ 
Eb = J2i eI^^ and e/ = We obtain Eqs.(|10|) and (|I2|) for the fields in the leading order, 

Hoipo = eo^o , (21) 

We first solve Eg. ([2l|) . Taking into account that 0o(±oo) —>■ ±v, we obtain the normal- 
izable solution, 

V'o(«^) = ^e-^°(-)x+ , (23) 
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with Uq{w) = G dw' (f)o(w') and the normahzation factor Af = dw e~'^^°^'^\ Here, 

X± are eigenspinors of o"2 satisfying a2X± = =tx± x±X± = 1- We obtain the energy 
eigenvalue, 

60 = ^ , (24) 

and take only eo positive, i.e., u = 1. Using Eg. ([23|) , we obtain the leading order fermi 
energy, 

Ef = Y. i^lHoi^^ = E I = R E K{2K) 

KM-^~°° KM ^ ^ K=l 

{Kmax + l){'2Kmax + I) 

= 3^ ■ ^''^ 

Here, Kmax is determined by the total fermion number, 

dw iPI^^ = J2 2K = KmaxiKmax + I) = Nf . (26) 
ft ill K=l 

From Eqs.(p5D and (p6|), we get 

(0) 2Arf / 1 \5 2nI Nf , ^ 

Ef = — L 1 + ~ — L + ^ for iV. > 1 . (27) 



The first term is the leading contribution to the fermi energy, which is the same as that 
obtained within the exact thin- wall approximation [0. The second term in Eq.(p7D is the 
correction caused by the effect of quantumizing the angular momentum, which we call AEf 
for the later discussion. 

We next consider Eq.(P^. Since the leading solution for the fermion satisfies ipla^ipQ = 0, 
the equation ( [2^ ) becomes 

2? = 5^o(^S-=) . (28) 

We know that the solution to the above equation satisfying (/)o(±oo) ±v is a kink, 

w 

<po{w) = t;tanh- , (29) 
db 

with 5b = 2/{-\f\v). Using Eq.(^), we obtain the leading order surface energy. 



with S = 2\/\v^ /?,. This also coincides with the result obtained in Ref.[^ within the thin- 
wall approximation. 



Combining Eq.(|30D with the first term in Eg. (pTl) , we get the leading order Fermi Ball 
energy, 

3 

Eo = ^ + AttR'J: . (31) 

( Note that we neglect the false vacuum volume energy, assuming that it is negligibly small. 
We discuss its magnitude in the cosmological context in Sec.|VT[) Minimizing Eg. (pi]) with 

respect to R, we get 



R = ^ , = ^ , \ , (32) 

(127rS)3 (87r)3A6w ^ 

which yields 

Eq = {87i)h^NfV . (33) 

Here, let us examine the stability of the Fermi Ball against the fragmentation using the 
leading order Fermi Ball energy. We compare two states: the one in which a single Fermi 
Ball has the fermion number Nf and the other in which n Fermi Balls have the fermion 
numbers less than Nf but keep the total fermion number to be Nf. Since the energy of the 
Fermi Ball in the leading order approximation is proportional to the total fermion number, 
the two states have the same energy. The leading order estimation cannot tell whether the 



Fermi Balls with large Nf produced in the early universe survive until now or not pj]. In 
order to examine the stability of the Fermi Ball, we calculate the higher order corrections 
in E in the next section. 

IV. HIGHER ORDER CORRECTION OF 5;,//?-PERTURBATION 

The next-to-leading order components of the fields satisfy 



Ho-^)^l = -{H^-e^)iJo = -Hl^Po , (34) 



K 

- A,3,^ _ ^ + ^ E (^J.W-. + W.sV-o) . (35) 

We here note 

/+00 1 r+oo 

dw Vo^i^o = 77 / we-^^o^'") = . (36) 

In order to solve Eg. (p^ , we first write ipi as the following form: 

i„{w) = {Uw)^^''\w) + , (37) 



where ip^"-\w) and ip^''\w) are the hnearly independent solutions to Eg . (|2T|) with the same 
eigenvalue eo, 



with W{w) = Jodw' e2^o("''). The equations (H) and {^) give 

— + — VTH— = , 
aw K aw 

^ + ^^e-^^o(-) = . 
dw 

The solutions to the above equations are 



i?3 Jo 

2K /•+°° 



du;' w'W{w')e-^^''^'"'^ -Ui{w) 



with Ui{w) = G Jq dw'(j)i{w'). We thus obtain the solution for ipi, 



where 



c+{w) 



c_ w 



2K' 
i?3 



fw poo 

/ du;'e2^° / dw" w;"e-2^°-f/iH 

JO Jw' 



K 



dw' w'e-^^o . 



The w-dependence of the component c± is shown in FIG.|l]. 



We next solve the Eq. (|35D . Replacing the fermionic source term of Eq. 
V-iVs^o = 6/(2Ar), we get 



dw AttRM 



KM 



h{w) 



(38) 
(39) 

(40) 
(41) 

(42) 
(43) 

(44) 

(45) 
(46) 



by ^/'(Jo-g^i 



(47) 



where 0i = The solution to the above equation is given by 



01 (W) 



cosh 



w 



dw' cosh^ 
a 



w 



dw' cosh — 



cosh^Y^ Jo Sb Jo cosh^^ 



Jo ) 
„ h{w") 



cosh^^^ 



dw" 



(48) 
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FIG. 1: The a2 = +1 (Figure (a)) and a2 = —1 (Figure (b)) component of the next-to-leading 
order solution for the fermion field with the Yukawa coupling constant G = 1. In the figures, c± 
and w are rescaled as ^JNfC± and w/5b, respectively, where Nf is the total fermion number of the 
Fermi Ball and 6^ is the domain wall thickness. The origin of the w axis denotes the center of the 
domain wall. Figure (a) shows that c+ is getting more centrally localized around the origin as the 
scalar self-coupling constant A decreases. Figure (b) shows that the width of c_ around the origin 
decreases with the decreasing A. 
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where we use 



°° , h(w) 

aw o — 

cosh^^ 



oc / aw — — n[w) = U 



/o cosn"^ Jo dw 

Ob 

Note that 0i(w) satisfies 0i(O) = and 0i(±oo) (see FIG.||). 
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(49) 
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FIG. 2: The next-to-leading order solution (pi for the domain wall field (p. In the figure, (pi and w 
are rescaled as \/^(pi = ^fXRcpi and w/5i,, respectively. Here, A is the scalar self-coupling constant 
and R and di, are the Fermi Ball radius and the wall thickness, respectively. The origin of the w 
axis denotes the center of the domain wall. The figure shows that (pi is not smooth at the origin 
in the limit of A — > 0, which coincides with the case where the spreading width of the fermionic 
source term in Eq.(R|) is neglected (see RehB). 



We note that the first order energy corrections E^p and El''' vanish since the integrands 
for them are odd functions of w. In the next-to-leading order, the fermi energy is 



dw 4H2^o + E / ^Hii^ 



2N] 



dw w^e-^^o + ^ 



hoo pw 
-2Uo 



dw we 



dw' (pi{w') 



3RWJ~oo 3RW 

5 

/ dw we-2^« / dw' e2^« / dw" w"e-^^^ for Nf > 1 

'Af J-oo Jo Jw' ^ 
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/ dw w^e H —7 / dw we / dw 0i(w') 

j-00 3A/ J-00 JO 



SAT 



and the surface energy is 



dw w;e-2^« / dw' e^^o / dw" w"e-^^'' , (50) 

-00 ^0 Jw' 



~oo 



2 




G 



2RN 

-2'k\^v^ / du; ^ 



/+00 — , 
dw 01^2 
-co 



cosh^l 



■ / du7 we / dw' 01 (w ) , (51) 

J-00 Jo 



SAT 

where R is replaced by Eq.(32). Taking into account Eq. (|32|) and AEf, the second term in 



the r.h.s. of Eg . (^71) , we obtain the Fermi Ball energy in the next-to-leading order, 

E2 = Ef +Ei^^ +AEf = CiX,G)v . (52) 

We see that C(A, G) does not depend on Nf] this is crucial to the discussion on the stability 
of the Fermi Ball below. Let us consider a single- Fermi Ball state and a n-Fermi Balls state 
both with the total fermion number taken to be Nf. The two states have the same leading- 
order energies Eq, while they have the different higher-order corrections E2: the former 
has a correction Cv and the latter nCv. Therefore, when C is positive, the former state 
has the lower energy and the Fermi Ball is stable against the fragmentation. We evaluate 
E2 numerically and find that C(A, G) is positive in a limited parameter region of A and G 
(see FIG.^. We see in FIG.|] that rather large value of G is allowed for the Fermi Ball 
to be stable. This situation, however, is much changed if the fermions have more degrees 
of freedom, e.g., belonging to a large multiplet of the internal symmetry and if the scalar 



11 



10 



E2=Q 

Stable region 

E2 > 



E2<0 
Unstable region 



1 1 — — — — I 

0.001 0.01 0.1 1 10 

X 



FIG. 3: The region of the scalar self-couphng constant A and the Yukawa coupUng constant G 
where the neutral Fermi Ball is stable (shadowed region) or unstable (blank region) against the 
fragmentation. The energy E2 in the figure denotes the deviation from the leading order energy 
of the Fermi Ball obtained within the thin-wall analysis. We see that rather large value of G is 
allowed for the Fermi Ball to be stable (see also FIG]^ and its caption). 

field (f) belongs to a singlet (see the forthcoming paper). Consider the fermion multiplet 
{1 < i < n) coupling to through the common Yukawa coupling constant G, and assume 
for simplicity the fermion number A^j = in common for each flavor i. In such a case, we 
see that the energy correction E2 is independent of but depends on n, and that a smaller 
value of G can stabilize the Fermi Ball (see FIG.^. Here, we emphasize that there is a 
certain region of the parameters where stable Fermi Balls are allowed to exist. 
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FIG. 4: The allowed region (shadowed) of the scalar self-coupling constant A and the Yukawa 
coupling constant G for the Fermi Ball with multi-fermions ^'j (1 < i < n) to be stable against the 
fragmentation (see text). The allowed region is illustrated for n = 1,3 and 10. The figure shows 
that the region extends as n is increases. 

V. COSMOLOGICAL OR OBSERVATIONAL CONSTRAINTS ON REGION OF 
PARAMETERS 

We very roughly examine here cosmological or observational constraints on the parameter 
region for the neutral Fermi ball. We define the parameter k = 27i^^^\^^^v which satisfies 



(53) 



with the Fermi ball mass Mf and the radius R being given by Eq.(^) and Eg. (P2|) , respec- 
tively. 

We first consider the energy density difference A which we have neglected so far. We 
know that A should be larger than the critical value A^ ~ /{lAA-K^Mpi) in order to avoid 
the domination of the black holes made up by domain walls in the total energy density of 
the early universe Such a finite A gives a volume energy = to the false 
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vacuum, having an effect to destabilize the Fermi Ball. As far as the volume energy is small 
enough, i.e., Ey = AnAMy^ /{Sk^^"^) < E2, the Fermi Ball is absolutely stable. This gives a 
constraint. 



,26 (GeVy fE,\'^ 



Mf < 2.9 X 10'^ [-^) GeV , (54) 

due to the condition for A to exist under the above constraints. 

We next consider the observational aspect of the neutral Fermi ball. If they are pro- 
duced in the early universe and have survived until present, they can contribute to the dark 
matter in the Galaxy. Let us assume here that the neutral stable Fermi ball has a sizable 
contribution to dark matter. We then have their flux E, 

Pdm uo in5/^GeV\ _2 -1 -1 
E ~ ~ 7.1 X 10 cm sec sr , (55) 

where pdm ~ 0.3 GeV cm~^ is the energy density of the dark matter in the Galaxy and 
Uo ~ 3 X 10^ cm sec~^ is the Virial velocity of the Fermi ball. We seek for the allowed region 
of the Fermi ball parameters through the use of currently available observational data. We 
use in the following the results of the experiments which searched for monopoles or heavy 
dark matter. 

In order to examine whether the Fermi ball can be detected in terrestrial experiments, 
we first consider the condition that the neutral Fermi ball should reach the detector passing 
through matter. The energy loss rate (energy loss per path length) is given by |jl2|, |13 

dE n 

— = p , (56) 

where a is the collision cross section with a nucleus in the medium, u is the velocity of the 
Fermi ball, and p is the density of the target matter. For simplicity, in the following we 
assume that the cross section is geometrically given by, 

a = 7TR^ . (57) 

From Eq.(^), the velocity u decreases exponentially with the path length, oc 
exp{—apx/Mf). Estimating the final velocity as Uc ~ 1.2 x 10^ cm sec""*^ |l^, we obtain the 
condition for the Fermi ball to reach the detector, 

Mf > O.UapL , (58) 
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where L is the path length. Substituting Eqs.(^) and (|57|) into Eq.(p8D yields 



pL 



1/3 



K > 4.6 X 10"^ — GeV . (59) 

\^gr cm~-^ j 

We next consider the efficiency of the detectors to observe the neutral Fermi ball. Let us 
examine the track detectors with mica and the scintillators in the MACRO fl^ and the 
KEK |jl9|] experiments. These experiments are sensitive to the neutral Fermi ball, if the en- 
ergy loss per path length is large enough, q{dE j pdx) > {dE/ pdx)min, where {dE/ pdx)min is 
the detection threshold for relativistic charged particles and q is the efficiency correction fac- 
tor. This condition with Eq.(^) gives cr > 1.7 x 10^^^[g~"'^(d-E'/pdx)mj„/GeV/gr cm~^] cm^, 
that is, 

q^'^ldE/pdx) 



M, = > 1.4 X 10 

TT 



GeV 



GeV/gr cm ^ 



GeV . (60) 



For the scintillator experiments in MACRO and KEK we take g = 1 for simplicity [|25 |. 
In case of the MACRO experiments, the conditions Eq. (|59D and Eq. (|60|) with pL = 3.7 x 
Iq5 q"^{dE/pdx)n,in ~ ^Imin ~ 6 McV/gr cm-2 give k > 3.3 GeV and 

Mf > 7.8 X 10^(K/GeV)3 GeV. The flux upper limit of F < 2.5 x 10"^^ cm-^sec^^sr-^ with 
Eq.(p5D gives the mass lower limit Mf > 2.8 x 10^^ GeV. In case of the KEK experiments, 
the conditions pL = 10^ gr cm~^ and q^^{dE/pdx) 

min ~ 0.01/mm ~ 20 KcV/gr Cm 

13 give K > 0.46 GeV and Mf > 2.6 x lO^^K/GeV)'^ GeV. The flux bound F < 3.2 x 



10-^^ cm-^sec^^sr^i gives Mf > 2.2 x 10^*^ GeV. 



Mf > 3.1 X 10^(K/GeV)3 GeV. The flux upper limit F < 2.3 x 10~20 cm-^sec-^sr-^ gives 
Mf > 3.1 X 10^5 GeV. 



For the truck detector with mica [|1^] we also take q = 1 for simplicity. In this case we 
take pL = 7.5 X 10^ gr cm~^ which is due to the fact that the mica was located at 3 km deep 
under the earth, and q^^{dE / pdx)min ~ 2.4 GeV/gr cm~^ They give n > 4.2 GeV and 

/ 
/ 

Recently very low background experiments are being done for the dark matter search 
|16| , ^ |2^ . Here, we analyze the results of CDMS experiments, which uses cryogenic Ge 
detector. The condition Eq.(|59D with pL = 2.6 x 10^ gr cm~^ gives k, > 1.2 GeV. The 
dark matter search experiments allow such a small cross section of the cold dark matter 
particle scattering with the target nucleus that the average number of collisions in the 
target is less than a unity {apd/rriGe ~ (ypd/Ace^p < 1, namely, a < 2.2 x 10^^^ cm^ 
with p = 5.3 gr/cm^, d = 1cm and Ace ~ 73). The experimental results are roughly 
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expressed as Fa oc cr/M/ < 3 x Ace x 10 '^^ cm^/GeV for Mf > 100 GeV, and this gives 
K > 3.9 X 10^ GeV. 

We illustrate the allowed region of Mf and k in Figure |^. In this figure the condition for 
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FIG. 5: The allowed regions (blank) of the Fermi ball mass Mf and the quantity k defined by 
K = 2tt^^'^\^/^v where A is the scalar self-coupling constant and v is the symmetry breaking scale. 
The upper shadowed region is excluded by the stability condition. The other shadowed regions are 



excluded by the experiments MACRO 0, KEK |T9|, MICA |T§ and CDMS [|T§. We take q = 1 
for simplicity (see text). 



the stability Eq.(^J) is also shown. One can see that there still remains the parameter region 
to be explored, especially the region of k ^ 10^ GeV for Mf ^ 10^^ GeV and k ^ 10^ GeV 
for 10^5 GeV < Mf < 10^^ GeV. 



VI. CONCLUSION AND DISCUSSION 

We have considered the neutral Fermi Balls in the thin-wall model where the domain wall 
thickness 6b is much smaller than the Fermi Ball radius R. In the case where the spreading 
thickness 6f of the fermion confined in the domain wall is negligibly small compared to 6b, 
the Fermi Ball is unstable against the fragmentation even if the finite 6b is taken into account 
0]. In the present paper, we have examined whether the Fermi Ball is stable or not if the 
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effect of the finite 5/ is included. 

In order to estimate the energy of the Fermi Ball, we have expanded the fields and the 
Hamiltonian in the power of R. At each level of the perturbation, we have examined the 
stability against the fragmentation. We have found that the energy correction in the next- 
to-leading order can stabilize the Fermi Ball in the limited region of the scalar self-coupling 
constant A and the Yukawa coupling constant G, as is shown in FIG.^. 

We have lastly given rough estimations for the allowed region of the parameters, Mj and 
K, for the neutral Fermi Ball to have a sizable contribution to the cold dark matter in case 
where the cross section of the Fermi Ball scattering with matter is of geometrical size. We 
have found that the allowed region is severely restricted by cosmological or observational 
constraints, but that there still remains the region open to the future exploration. 

In the present paper, we have dealt with the Fermi Ball in a semi-classical manner. 
The quantum corrections, such as the radiative correction coming from the scattering of 
the fermions in the domain wall, may affect the stability of the Fermi Ball. They will be 
discussed elsewhere. 
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